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ABSTRACT 

We present X-band data for the brightest cluster galaxies (BCGs) from the ESO Dis- 
tant Cluster Su rvey (EDisCS). These data are combined with the photometry published by 
lAragon-Salam anca et al. (1998) an d a lo w-redshift comparison sample built from the BCG 
catalogue of rvon der Linden et al] (|2007[) . BCG luminosities are measured inside a metric 
circular aperture with 37kpc diameter. In agreement with previous studies, we find that the 
i^-band Hubble diagram for BCGs exhibits very low scatter {^ 0.35 mag) over a redshift 
range of < z < 1. The colour and rest-frame iiT-band luminosity evolution of the BCGs 
are in good agreement with population synthesis models of stellar populations which formed 
at z > 2 and evolved passively thereafter. In contrast with some previous studies, we do not 
detect any significant change in the stellar mass of the BCG since z ^ \. These results do 
not seem to depend on the velocity dispersion of the parent cluster. We also find that there is 
a correlation between the ID velocity dispersion of the clusters (ctci) and the /-C-band lumi- 
nosity of the BCGs (after correcting for passive evolution). The clusters with large velocity 
dispersions, and therefore masses, tend to have brighter BCGs, i.e., BCGs with larger stel- 
lar masses. This dependency, although significant, is relatively weak: the stellar mass of the 
BCGs changes only by ^ 70% over a two-order-of-magnitude range in cluster mass. Fur- 
thermore, this dependency doesn't change significantly with redshift. We have compared our 
ob servational results w i th the hierarchical galaxy formation and evolution model predictions 
of iDe Lucia & BlaizotI (l2007b . We find that the models predict colours which are in reason- 
able agreement with the observations because the growth in stellar mass is dominated by the 
accretion of old stars. However, the stellar mass in the model BCGs grows by a factor of 3- 
4 since z — 1, a growth rate which seems to be ruled out by the observations. The models 
predict a dependency between the BCG's stellar mass and the velocity dispersion (mass) of 
the parent cluster in the same sense as the data, but the dependency is significantly stronger 
than observed. However, one major difficulty in this comparison is that we have measured 
magnitudes inside a fixed metric aperture while the models compute total luminosities. 

Key words: galaxies: clusters — galaxies: formation — galaxies: evolution — galaxies: 
elliptical and lenticular, cD — infrared: galaxies 
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1 INTRODUCTION 

Brightest cluster galaxies (BCGs) form a very special class of 
galaxy. They lie in local minima of the cluster potential well and 
are extremely massive and luminous, forming their own luminos- 
ity fu nction (Sandage 1976; Dressier 1978; Bernstein & Bhavsar; 
uOOlh. BCGs have been the focus of many studies (see , 
e.g.. iHumason et alj 1 19561; ISandagd 1 19721; lOunn & Okd |l974 



iHoessel et al.lll98(]| ; [oegerle & Hoessellll99ll ; lNelson et aljfzoa 
and are well known for the low scatter in their absolute magnitudes 
(~ 0.3 mag; see, e.g., Sandage 1988). Early work involving BCGs 
took advantage of this small scatter in absolute magnitudes, com- 
bined with the extreme luminosities of these galaxies, to examine 
the cosmological bulk flow out to higher redshifts and to try to de- 
termine the value of the deceleration parameter go using the BCGs 
as standard candle s (e.g. Sandage 1972; Sandage & Hardv 1973; 
iGunn & Okelll975l ; iLauer & Postmanlll992h . More recent studies 
of BCGs have focused on unravelling the formation and evolution 
of these objects in order to place constraints on galaxy formation 
models. Aragon-Salamanca, Baugh & Kauffmann (1998; hereafter 
ABK) examined the A'-band Hubble diagram for BCGs up to a 
redshift of ^ 1. They found that the BCGs did not exhibit any 
luminosity evolution over this redshift range, suggesting that the 
stellar mass of the BCGs had grown by a factor 2-4 since z ^ 1, 
depending on the co smology assumed. Collins & MannI ( 119981) and 
ISrough et alj (2003) also studied the AT-band Hubble diagram for 
the BCGs in an X-ray selected cluster sample. These authors found 
that the amount of stellar mass growth shown by the BCGs was 
dependent on the X-ray luminosity of the host cluster, with BCGs 
in high X-ray luminosity clusters showing no mass accretion since 
z ~ 1 and BCGs in low X-ray luminosity clusters growing by a 
factor of ~ 4 in thi s time. 

More recentlv lDe Lucia & Blaizoj 420070 have used a combi- 
nation of N-body simulations and semi-analytic techniques to study 
the formation and evolution of BCGs. They find that the formation 
history of local BCGs is extremely hierarchical, with half the mass 
of a typical BCG being locked up in a single galaxy after z ~ 0.5. 
The stars that make up the BCGs are formed very early in separate, 
small galaxies which then assimilate over time to form the BCG. 

In this paper we extend the work carried out by ABK, enlarg- 
ing their high-redshift sam ple with the BCGs from the ESO Distant 
Cluster Survey (EDisCS; Iwhite et al.|[20o3) . The local compari- 
son sample is al so improved usi ng the BCGs from an SDSS-based 
cluster sample jvon der Linden et al...2007.) . Moreover, additional 
information such as the velocity dispersions of the host clusters 
is now included in the analysis. In section 2 we present the data 
and describe our main observational results concerning the stellar 
mass growth of the BCGs, the evolution of their colours and the 
relationship between the velocity dispersion of the host cluster and 
the BCG luminosity. Section 3 compares the observational results 
with t he semi-analytic model predictions of i De Lucia & Blaizoll 
( 120071) . Section 4 summarises our main conclusion s. We assume 
Ho = 71kms"^Mpc~\ ilm = 0.27, Oa = 0.73 i lSpergel et al.] 
120030 throughout the paper. 



2 OBSERVATIONAL RESULTS 

2.1 The sample 

One part of the data analysed here is taken from th e ESO Dis- 
tant C luster Survey (EDisCS), presented in detail in IWhite et al.l 
( 120051) . A brief description is given here for completeness. The 



candidate clusters for EDisCS were chosen from among the high- 
est surface brightnes s objects in the Las Ca mpanas Distant Cluster 
Survey (LCDCS) of iGonzalez et alj bOOlh . Potential cluster can- 
didates were identified in two redshift bins, mid-z clusters (0.45 < 
2cst < 0.55) and high-z clusters (0.75 < Zcst < 0.85), where 
Zest is the estimate d cluster redshift base d on the magnitude of the 
putative BCG (see lGonzalez et alJuOOll for details). The potential 
clusters were checked carefully to ensure the detection appeared re- 
liable and deep optical and near-infrared photometry was taken for 
20 fields with confirmed cluster candidates. The optical photom- 
etry was taken using FORS2 at the VLT in direct imaging mode. 
The mid-z clusters were imaged in the B, V and / bands, the high- 
z clusters in the V, R and / bands. Typical integration times for the 
optical data were 45 minutes for the mid-z clusters and 120 minutes 
for the high-z clusters. The near-infrared J and K photometry was 
obtained using SOFI at the NTT (Aragon-Salamamca et al. 2008, 
in preparation). Two fields were not imaged due to bad weather and 
a further field was rejected as its spectroscopic redshift histogram 
suggested the cluster was more likely to be a chance superposition 
of objects spread out over a range of redshift. However, confirma- 
tion spectroscopy with FORS2 at the VLT showed that several of 
the confirmed fields were found to contain multiple clusters at dif- 
ferent redshifts, which compensated for the loss of the other 3 fields 
dGonzalez et al.l2002l ; rwhite et al.l200i) . The EDisCS sample con- 
sidered here contains K-band data for a total of 21 BCGs spread 
over a redshift range of approximately 0.4 to 0.96. (see Table [TJ. 
A minimum of 150 minutes of integration in the Ks band was ob- 
tained for the mid-z clusters and a minimum of 300 minutes for the 
high-z ones. HST/ACS data was also obtained for the 10 highest 
redshift clusters. Typical random errors for the ground based pho- 
tometry are ~ 0.007 mag. Errors in the photometric zero-points 
are ~ 0.02 mag. The EDisCS BCGs were identified by examin- 
ing each cluster individually and choosing the most suitable galaxy 
in each case, taking into account redshift, brightness and location 
of surface density peak. Further information re garding the EDisC S 
photometry and BCG selection is presented in lWhite et al.l ( 120051) . 

In addition to th e EDisCS BCGs we make u se of a sample of 
local BCGs studied in lvon der Linden et al.l J2007r) . Only a brief de- 
scription of how these BCGs were selected is included here. For a 
more detailed account, see von der Linden et al. (2007 ). The start - 
ing point of this sample was the C4 catalogue ( Mille r et al.l 20050 
of clusters in the Sloan Digital Sky Survey (SDSS; iStoughton et al.l 
1 20021 lYork et al.l [200(1). However, due to fiber collisions in the 
SDSS (no two fibers can be placed within 55" of each other), about 
a third of the BCGs are missed by the C4 algorithm, which is based 
on the spectroscopic data alone. A semi-automatic algorithm to 
identify the brightest, central galaxy o f these clusters is presented in 
detail in lvon der Linden et alj ( 120070 . Furthermore, an iterative al- 
gorithm is used to determine the redshift z, velocity dispersion CTci, 
and virial radius i?200 from cluster members within I-R200 from 
the BCG and within ±3crci of the cluster redshift. The final clus- 
ter sample contains 625 groups and clusters. Of their BCGs, 604 
are found in the 2M ASS Extended Source Catalogue (2MASSX; 
ISkrutskie et alj|2006r) . Further cuts to the local sample during the 
process of obtaining reliable A"-band aperture magnitudes (see Sec- 
tionim reduced the total number of SDSS/2MASS BCGs to 598, 
covering a redshift range 0.02 ^ z ^ 0.10. These galaxies will be 
used in the following sections to build a low-redshift comparison 
sample matching, as closely as possible, our high-redshift one. 

We also combine the EDisCS and SDSS/2MASS data with 
the photometry published by ABK. We refer the reader to the orig- 
inal paper for a description of their data. ABK provides A-band 
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Table 1. Corrected rest-frame X-band magnitudes for the EDisCS BCGs 



Cluster 


BCGID 




2c; 


Ko 


(T(,i(km/s)'* 


cll301-1139a 


EDCSNJ1301351- 


-1138356 


0.3969 


16.0 


391±ii 


C11202-1224 


EDCSNJ1202433- 


-1224301 


0.4240 


16.7 


518±?g, 


cll037-1243a 


EDCSNJ1037523- 


-1244490 


0.4252 


16.9 


537±ti 


C11059-1253 


EDCSNJ1059071- 


-1253153 


0.4564 


16.0 


510±ii 


C11018-1211 


EDCSNJ1018467- 


-1211527 


0.4734 


16.3 


486±ii 


ell 138-1 133 


EDCSNJ1138102- 


-1133379 


0.4796 


17.0 


732±|i 


C11301-1139 


EDCSNJ1301402- 


-1139229 


0.4828 


16.4 


687±ii 


C11420-1236 


EDCSNJ1420201- 


-1236297 


0.4962 


16.1 


218±tg 


C11411-1148 


EDCSNJ1411047- 


-1148287 


0.5195 


16.2 


7io±;ii 


C11232-1250 


EDCSNJ1232303- 


-1250364 


0.5414 


16.3 


w%otll'' 


ell 037- 1243 


EDCSNJ1037514- 


-1243266 


0.5783 


16.4 


319±ii 


C11353-1137 


EDCSNJ1353017- 


-1137285 


0.5882 


16.5 


(>(>6t\ll 


clll03-1245a 


EDCSNJ1103349- 


-1246462 


0.6261 


17.4 


336±lg 


C11227-1138 


EDCSNJ1227589- 


-1135135 


0.6357 


\1A 


574±|i 


C11054-1146 


EDCSNJ1054244- 


-1146194 


0.6972 


17.2 


589±?e 


clll03-1245b 


EDCSNJ1103365- 


-1244223 


0.7031 


17.1 


252±ii 


C11040-1155 


EDCSNJ1040403- 


-1156042 


0.7043 


17.5 


418±ti 


C11054-1245 


EDCSNJ1054435- 


-1245519 


0.7498 


17.5 


504±ii3 


C11354-1230 


EDCSNJ1354098- 


-1231015 


0.7620 


17.4 


648±1?§ 


C11216-1201 


EDCSNJ1216453- 


-1201176 


0.7943 


17.0 


1018+11 


C11103-1245 


EDCSNJ 11 03434- 


-1245341 


0.9586 


18.0 


534±;oj 




HHallidav et alJi2004l),lMilvan2- 


Jensen etalJbOOSi) 
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Figure 1. Top: histogram of the cluster velocity dispersions for the complete 
SDSS/2MASS (solid line) and the high-z cluster samples (dashed lines). 
Bottom: as above, but using one of the SDSS/2MASS tT^i-matched sub- 
samples. The high-2 cluster histograms have been normalised to have the 
same number of objects as in the SDSS/2MASS samples. 



photometry for a sample of 25 BCGs covering a redshift range of 
0.023 < z < 0.92. These galaxies will also form part of our study. 



2.2 The K-band Hubble diagram 

In order to allow a direct comparison between the data sets, we 
closely followed the method adopted by ABK when construct- 
ing the K-b&nA Hubble diagram, changing only the assumed cos- 
mological parameters. The magnitudes of the BCGs were mea- 
sured including all the light contained inside a fixed metric circular 
apertu re centred on the BCG itself (cf. lSchneider. Gunn & Hoessell 
119831) . ABK used a metric aperture of 50kpc diameter, assuming 
Ho — 50kms^^ Mpc^^, and carried out their analysis for both 
go = 0.0 and f/o = 0.5. We have assumed the cosmological pa- 
rameters from lSpergel et alj ( 120031) as derived from the Wilkinson 
Microwave Anisotropy Probe (WMAP), Hq = 71kms~^ Mpc"^, 
Sim = 0.27, Q.K = 0.73. An aperture of 50kpc in the go = 0.0 
cosmology assumed by ABK is equivalent to an aperture of diame- 
ter ~ 37 kpc in the WMAP cosmology for the redshifts considered 
here. 

For the EDisCS BCGs, the aperture magnitudes measured di- 
rectly from the i mages were corrected for galact i c extin ction using 
the dust maps of ISchlegel. Finkbeiner & Davij ( Il998l) . In the K- 
band these corrections were always small (typically ~ 0.01 mag) 
since the galactic latitude of the fields is relatively high. A k- 
correction was also applied following the sa me method used b y 
ABK using the fc-correction published by Man nucci et al.l ( 120011) . 
The combined uncertainties in the derived K-hanA magnitudes are 
~ 0.05-0.08, dominated by the uncertainties in the fc-correction 
(cf. ABK). 

The method employed to measure the required aperture mag- 
nitudes for the SDSS/2MASS BCGs is as follows. We searched 
for the positions of the 625 objects selected from the SDSS BCG 
sample in the 2M A SS Extended Source Catalogue (2MASSX; 
ISkrutskie et alj|2006h . Twenty-one of the objects were not found, 
reducing the SDSS BCG comparison sample to 604 objects. 
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2MASSX provided magnitudes taken in circular apertures of in- 
creasing radius (5", 7", 10", 15", 20", 25", 30", 40", 50", 60" 
and 70"). The maximum aperture for which a magnitude was avail- 
able varied from object to object. These aperture magnitudes were 
then interpolated to estimate the value corresponding to a physical 
aperture of 37 kpc at the redshift of each BCG. Unfortunately, for 6 
of the SDSS BCGs the required aperture was larger than the maxi- 
mum aperture that was available from the 2MASSX database. For 
these objects it would have been necessary to extrapolate to esti- 
mate the required aperture magnitudes. Rather than using uncertain 
extrapolated magnitudes, we decided to discard these objects, leav- 
ing a final sample of 598 SDSS/2MASS BCGs. Simple statistical 
tests show that these galaxies represent a random subsample of the 
original SDSS/2MASS BCG sample, unbiased with respect to their 
redshift distribution, velocity dispersion or luminosity. Galactic ex- 
tinction corrections were applied to these objects in the same way 
as for the EDisCS BCGs. The correction was in general small (typ- 
ically ~ 0.01 mag). A fc-correction was also applied to the data as 
before. 

Note that both the 2MASS and the EDisCS photometric data 
are in the K^ system, while the ABK data are in the UKIRT K 
system. However, the differences between the photometric systems 
are ^0.01 mag and thus negligiblqj. 

In section l2!4l we will see that the luminosity of the BCG de- 
pends significantly on the velocity dispersion (ctci) of the parent 
cluster. Thus, in order to make meaningful comparisons between 
the low-redshift and high-redshift BCGs, the velocity dispersions 
of the cluster samples need to be reasonably well matched. The top 
panel of figure [T] shows that the SDSS/2MASS sample contains a 
much larger proportion of low ad clusters than the high-z sam- 
ple, thus a direct comparison cannot be made. However, given the 
number of clusters in the SDSS/2MASS sample, it is possible to 
build reasonably-large subsamples of low-2 clusters with CTcI dis- 
tributions that are well matched to that of the high-z sample. A 
set of 1000 local cluster subsamples, each with a CTcI distribution 
matched to our high-z sample, was built by randomly-selecting 
SDSS/2MASS clusters in the right proportion. Each of the CTcI- 
matched subsamples contains 50 SDSS BCGs. A typical realisation 
is shown in the bottom panel of figure [T] In this sample-matching 
process we did not take into account the redshift evolution of Cci 
because such approach would have been model d ependent. How- 
ever, the expected evolution is small enough (e.g.. |Poggianti et al.l 
l2006h and the dependence of the BCG's luminosity on Uci weak 
enough (cf. section l2!4t that our conclusions would not change if 
we had taken such evolution into account. 

In the analysis described below, we carried out comparisons 
between the high-z sample and the 1000 realisations of the low-z 
(Tci-matched sample. The quantitative results we quote always refer 
to the mean of all these comparisons. For display purposes, and for 
the sake of clarity, we sometimes only plot the points corresponding 
to one of these realisations whose properties are close to the mean. 

Figure |2] shows the rest-frame A'-band Hubble diagram for 
the BCGs in the SDSS/2MASS, EDisCS and ABK clusters. The 
estimated random uncertainties in the A'-band magnitudes add little 
to the scatter of the points, as they are ~ 0.05-0.08 mag for the 
EDisCS and ABK BCGs and ~ 0.09 mag for the SDSS/2MASS 
BCGs. 

The solid line in figure |2] is the no-evolution prediction. 



1 See the Supplement to the 2MASS All Sky Data Release, 
|http://www.ipac.caltech!edu/2mass/releases/allsky| 



Table 2. Corrected rest-frame K-hand magnitudes for a typical a^i- 
matched SDSS subsample. The properties of this subsample are close to 
the mean of the 1000 subsamples (See section |2!4l ) 



SDSS BCG ID 


2cl 


Ko 


o-cl(km/s)'' 


587726014532550731 


0.02015 


8.7 


502+ g 


587722953306931328 


0.04362 


12.0 


264±ii 


588015507664928793 


0.04437 


10.5 


933±?i 


587728932421042326 


0.04631 


11.2 


556±il 


587728670418272398 


0.04845 


11.5 


570±g 


588848899917611022 


0.04876 


10.9 


mtil 


587729387149852828 


0.05189 


11.7 


547±fgo 


587727227305263139 


0.05558 


10.7 


903±ii 


587722952231944566 


0.05566 


12.9 


685±ii 


587734304344178742 


0.05625 


12.4 


9151?^ 


587733081347195107 


0.05643 


11.4 


441111 


587722983377535296 


0.06003 


11.9 


506±tf 


587729751131029774 


0.06029 


11.2 


512±i? 


587732134304350396 


0.06288 


11.7 


613j;S 


587731185668849907 


0.06517 


11.9 


io84±;;i 


588010360151146773 


0.06901 


12.2 


motii 


587733081882755176 


0.07107 


13.0 


543±i| 


587726032259973170 


0.07174 


11.2 


sutii 


587732582591168714 


0.07267 


12.1 


582±18° 


587735666921898155 


0.07371 


11.4 


670±™ 


587731680111820927 


0.07587 


12.5 


902+?67 


587727225689538702 


0.07618 


11.5 


812±tf 


587726016681476254 


0.07738 


11.9 


966±i§ 


587727212807192844 


0.07868 


12.8 


765±?L 


587726031184920689 


0.07919 


13.1 


327±?1 


587725576423342431 


0.08012 


12.0 


1156±iH 


587726031184330907 


0.08044 


13.0 


677±ii 


588015508205994134 


0.08157 


12.0 


497±ii 


587725505559986486 


0.08182 


12.6 


613±i^ 


587730817361641713 


0.08255 


13.1 


406+ ii 


587726877271130340 


0.08257 


13.4 


315±il 


587726032776593461 


0.08316 


12.6 


424±™ 


587733080814780480 


0.08330 


12.3 


625±?1 


587730816824508498 


0.08512 


12.5 


22911^4 


587729782810345782 


0.08578 


12.3 


530±j;i 


588013381670994051 


0.08581 


12.6 


mtii 


587729752747278348 


0.08803 


11.7 


951±S 


587733410986590314 


0.08804 


11.9 


105tl% 


588848899917873272 


0.08820 


12.2 


1020+iI 


588848900454678767 


0.08856 


11.9 


708±JSi 


587733604808458243 


0.08901 


13.1 


556±?;o 


587736618787733576 


0.09010 


12.4 


974±if 


587722982834438274 


0.09059 


12.2 


V48±ii 


587726100415774941 


0.09243 


12.3 


513±ii 


587735694836039706 


0.09356 


13.3 


554±;?= 


587727212274057247 


0.09465 


12.4 


513±fg5 


588010359621156952 


0.09742 


13.2 


1074±}°i 


587726032234741909 


0.09767 


12.5 


388±;;i 


587729407545639093 


0.09903 


12.1 


until 


587726032245620919 


0.09947 


13.3 


746±iS 
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which takes into account dimming solely due to distance modu- 
lus as a function of redshift. The other two lines plotted are pas- 
si ve evolution predic t ions f rom the population synth esis models 
of iBruzual & Charloll ( l2003h , using a lSalpeteil h955l) initial mass 
function with solar metallicity. We have assumed a fixed metallic- 
ity for the BCGs as the actual metal content and its evolution is 
unknown. However, ABK showed that changing the metal content 
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Figure 2. Top plot: the apparent rest-frame ii'-band magnitudes versus redshift for the BCGs. The EDisCS points are shown as plus signs, the data from 
ABK are represented by circles, and the galaxies in one of the SDSS/2MASS Cci-matched subsamples is shown as crosses. This SDSS/2MASS subsample 
has properties close to the mean of the 1000 realisations. Small dots coiTespond to the complete SDSS/2MASS BCG sample, shown here for comparison 
(see text for details). The magnitudes have been measured in a fixed metric aperture of 37 kpc. The solid fine shows the no-evolution prediction, and the 
dashed and dotted lines show the predicted evolution for passively-evolving stellar populations formed at redshifts of 2 and 5 respectively. The model lines 
are normalised at low-z to the mean of the 1000 CTd-matched SDSS/2MASS subsamples. Bottom plot: the top panel shows the data with the no-evolution 
prediction subtracted (symbols as above). The middle and bottom panels show, respectively, the data after subtracting the passive evolution predictions for 
•^form = 2 and Zform = 5. The dashed line in each panel indicates the average least-squares fits to each of the 1000 low-z subsamples plus the EDisCS and 
ABK BCGs. 



makes little difference to the passive luminosity evolutiorQ The 
dashed line in figure [2] shows the predicted evolution of a stellar 
population which forms at a redshift of 2 and evolves passively 
thereafter. Similarly, the dotted line is the prediction for a passively 
evolving stellar system formed at a redshift of 5. The prediction 
lines are normalised at low-z to the mean of the 1000 crd-matched 
SDSS/2MASS subsamples. Contrary to the findings of ABK, who 
found no evidence for luminosity evolution in the if-band lumi- 
nosity of BCGs, figure |2] indicates relatively small, but significant, 
luminosity evolution: BCGs at z ~ 1 appear to be, on average, 
1.7 times brighter in the rest-frame K band than z ~ ones. In- 



^ The metallicity of the BCGs is discussed in section |231 



deed, the BCGs follow quite closely the passive luminosity evolu- 
tion lines, with a scatter of ~ 0.35 mag. The different behaviour 
we have found is mainly due to the low-z: comparison sample used 
to derive the low-z normalisation of the Hubble diagram. Due to 
the limited amount of K-band data for BCGs available at the time, 
ABK used a small ad-hoc sample of only 4 low-z BCGs. The large 
uncertainty on the normalisation from such a small sample was 
exacerbated by the fact that most of these clusters had relatively 
high velocity dispersions, biasing their magnitudes. This highlights 
the dangers of comparing low- and high-z cluster samples without 
carefully matching their properties. Indeed, had we used the entire 
SDSS/2MASS BCGs sample (small dots in Figure^ to normalise 
the Hubble diagram, we would have derived a much stronger evo- 
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lution since the majority of the local clusters have low velocity dis- 
persions and their BCGs are therefore too faint (cf. section l2!4] l. 

The three lower panels of figure |2] show the same data after 
subtracting the no-evolution and passive-evolution predictions. As 
mentioned above, it is clear from the top panel that the BCGs show 
a relatively small, but significant, luminosity evolution over the red- 
shift range plotted. The results of subtracting the passive evolution 
model predictions from the data are shown in the middle and bot- 
tom panels assuming ztorm = 2 and Zform = 5 respectively. The 
lines of best fit in these panels (derived from the average of the 
1000 low-z subsamples) are very close to the passive luminosity 
evolution predictions themselves. This is in agreement with the ob- 
served colour evolution of these objects, which suggests that their 
stars were formed at a redshift z > 2 and then evolved passively 
(see section |23] >. 

If the stellar populations of these objects are passively evolv- 
ing, as their colour evolution suggests (section [23] ). the observed 
luminosity evolution can be used to place limits to the amount 
of growth in the stellar mass of the BCGs. Following the method 
of ABK we parameterize the stellar mass evolution as Af*(z) = 
Mi,{Q) X (1 + 2)^. Least-squares fits to the data points in the bot- 
tom panels of figure[2]give the following results: 



2form — O 



7 = -0.3 ±0.2 
7 = +0.1 ±0.2 



This means that we do not detect any significant change in the stel- 
lar mass of these objects. Our data is compatible with zero growth. 
Formally, the derived 7 values and their associated errors imply 
that, if the stellar populations in these galaxies formed at Zform = 2, 
the stellar mass cannot have grown by more than a factor 1.9 be- 
tween z = 1 and 2 = (3-(t upper limit). For Zform = 5, the upper 
limit reduces to a factor of 1.4 (i.e., only 40%). This contrasts with 
the results obtained by ABK, who suggested stellar mass growths 
for the BCGs of factors between 2 and 4 since 2 ~ 1. As men- 
tioned above, this difference is mainly due to their low-2 normali- 
sation. We compare the observed mass growth with the results from 
sem i-analytic models i n sect ion[3] 

ICollins & Mmr] j 19981) and iBrough et al] ( l2002h found that 
the amount of stellar mass growth shown by the BCGs could be 
dependent on the X-ray luminosity of the host cluster, with BCGs 
in high X-ray luminosity clusters showing little or no stellar mass 
growth since 2 ~ 1, while BCGs in low X-ray luminosity clusters 
exhibit a large variation in their stellar mass growth, from clusters 
showing no growth to a few growing by factors of up to ~ 4. As 
discussed in i ]2.4| the cluster X-ray luminosity and velocity disper- 
sion (Jci can be taken as a proxies for cluster mass. If we divide 
our low- and high-2 cluster samples into two equal subsamples, 
one with a^i ^ 654kms~^ and the other with CTci > 654kms~^ 
we find that the inferred stellar mass growth does not seem to de- 
pend on (Tci: the stellar mass of the BCGs seems to have remained 
roughly constant for both subsamples (albeit, with weaker formal 
limits given the smaller sample sizes). This is shown in figure [3] 
which presents again the lower panels of figure[2]with the data split 
into high and low a^i subsamples. There is no significant difference 
in the luminosity evolution rate of the two ad subsamples. How- 
ever, there is clearly an offset between the lines of best fit in each 
panel of figure [S] with the high-crd sample being brighter in each 
case. This implies a significant dependence of the luminosity of the 
BCG on cluster velocity dispersion, as discussed in section 2.4. 

To test whether our conclusions depend strongly on the choice 
of galaxy evolution models, we have re-calculated the expected 
passive A"-band luminosity evolution (cf. figure[2ll using the models 




log(l+z) 

Figure 3. The data from the K-haaA Hubble diagram with the model pre- 
diction lines subtracted away (cf. bottom panels of figure|2). The data has 
been divided into high- and Iow-cTj-i subsamples. The high-(T(.i subsample 
is plotted as solid square points and the low-o-d subsample as open circles. 
The dashed and dotted lines in each panel represent least-squares fits to the 
high-cTci and low-cr^i data respectively. See text for details. 



of iMaraston et alj yOOg). Since these models use a very different 
prescription for TP-AGB stars, the near-IR luminosity evolutio n is 
significantly different from that of Bruzual & CharlotI (120030 for 
stellar populations with ages ~ 0.5-1 Gyr. However, for the range 
of ages relevant to our analysis (2.6-12.5 Gyr), the model predic- 
tions agree within ~ 0.1 mag for a given metallicity and IMF. We 
are thus confident that the conclusions of this section are reasonably 
model-independent. 



2.3 Colours 

Figure|4]shows the observed-frame V — I and I — K colours for the 
EDisCS and SDSS/2MASS BCGs, plotted as a function of redshift. 
For the EDisCS BCGs a 2" radius circular aperture was used to 
measure the colours. This aperture is large enough to sample a sig- 
nificant fraction of the BCG's light, while keeping sky-subtraction 
errors smalo. The 2" radius at the average redshift of the EDisCS 
BCGs, (~ 0.6), corresponds to an aperture of radius ~ 10" at the 
average redshift of the SDSS BCGs (~ 0.07). Optical and near- 
infrared photometry for the SDSS/2MASS sample was measured 
inside circular aperture s with 10" radius, and tra nsformed into the 
VIK system following Isianton & RoweisI (^2007). 

The colours were corrected for galactic extinction but not k- 
corrected. The lines shown in figu re [4] are predictions create d us- 
ing the population synthesis code of lBruzual & CharlotI ( 120031) . The 



^ Measuring galaxy colours inside the metiic aperture used in section 12.21 
(average radius^^ 2.5") makes no significant difference to the conclusions 
of this section. 
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solid lines represent the expected col our evolut i on of a stellar pop- 
ulation with solar metallicity and a ISalpeteii ( 119550 IMF which 
formed in an instantaneous burst at 2 = 2 and then evolved pas- 
sively. The dotted and dashed lines show the expected evolution 
for similar stellar populations that formed at 2 = 5 and z = 1 
respectively. It is clear from figure |4] that, within the model uncer- 
tainties (see below), there is reasonably good agreement between 
the data and the passive evolution models with 2;form = 2 and 5. 
As expected, the model for Zform = 1 does not fit the data at all. 
The agreement between the colours of the BCGs and the passive 
evolution models is a strong indication that the BCGs are com- 
posed of stellar populations that formed at z > 2 and have been 
evolving passively since. Similar conclusions are obtained using 
the .Maraston et al.. (2006.) models. For the range of ages relevant to 
2form betwee n 2 and 5, these models pr edict colours which agree 
with those of i Bruzual & CharlotI ( l2003h within ~ 0.1 mag. This 
should be taken as an indication of the typical uncertainties asso- 
ciated with model colours. For that reason, the fact that the V — I 
BCG colours seem to agree with the Zform = 5 models while for 
I — K there is better agreement with Zform = 2 is not significant. 
Nevertheless, since the V band samples the rest-frame ultraviolet, 
a very small amount of recent star formation (say, ~ 1% by mass) 
could easily make the observed V — I colours bluer by ~ 0.1 while 
having negligible effect on I — K. 

As the metallicity of the BCGs is unknown, the colour evolu- 
tion predictions were also calculated for a stellar population with 
approximately double solar metallicity and also approximately half 
solar metallicity. These predictions were found to lie nowhere near 
the plotted data points, being far too red and blue respectively 
to reproduce the colour evolution exhibited by the BCGs. This 
shows that although the colour evolution predictions are sensitive 
to changes in metallicity, the assumption of approximately solar 
metallicity is a reasonable one. 

2.4 BCG A' -band luminosity and cluster velocity dispersion. 

Arguably, the most fundamental property of a cluster is its mass, 
so it is not unreasonable to expect that the luminosity (and stel- 
lar mass) of the BCGs may be related to it. Several observational 
properties can be used as proxies for the cluster mass, including the 
cluster X-ray luminosi ty and temperature and the cluster velocity 
dispersion ad- Indeed. lEdge & Stewarli ( Il99lh found a significant 
correlation, in a low-« sample, between the optical luminosity of 
the BCGs and the X-ray luminosity and temperature of the parent 
clusters in the sense that hotter and more X-ray bright clusters host 
brighter BCGs. 

Velocity dispersions are available for all the EDisCS and 
SDSS clusters and 10 of the clusters studied by ABK. Columns 4 of 
tables[T]and[2]show the velocity dispersions CTci of the EDisCS and 
SDSS/2MASS (Tci-matched cluster samples. A full description of 
th e calculation of the ctc I value s for the SDSS clusters can be found 
in I von der Linden et alj (12007). Information regarding th e EDisCS 
(Tci va lues can be found in Halli dav et al. (2004), Pogg ianti et al.l 
( I2OO 6I and Mi lvang- Jensen et alj (12008). The veloci ty dispers ions 
for the ABK cluste rs are taken from lGirardi & Mezzettil 1 1200 ih and 
iLubin et al.l j2002i) . Note that in all these studies the cluster velocity 
dispersions have been determined following very similar methods, 
and are thus directly comparable. 

In order to test for any relationship between the BCG A'-band 
luminosity (or stellar mass) and the cluster mass, we have calcu- 
lated the residuals around the three lines of best fit shown in the 
bottom three panels of figure|2]and plotted them against logj^g (uci ) 
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Figure 4. Plot of V — I (crosses) and I ~ K (dots) observer-frame 
colours for the EDisCS BCGs (2 > 0.4) and the SDSS/2MASS BCGs 
(z < 0.4). The lines are p redictions from the population synthesis code 
of lBruzual & Charloj J20O3l) . The solid lines represent the colour evolution 
prediction for a stellar population formed at a redshift of 2. The dotted and 
dashed lines are the predictions for stellar populations with Zform = 5 and 
•^form = 1 respectively. The predictions were computed using a Salpetea 
119551) IMF and solar metallicity. 



in figure [7] This procedure should correct for redshift-related ef- 
fects. The panels show, from top to bottom, the residuals around 
the line of best fit after the rest-frame A'-band magnitudes have 
been corrected for no evolution, passive evolution with Zform = 2 
and passive evolution with Zform = 5, respectively (cf. figure |2]l. 
The symbols for the EDisCS and the ABK samples are the same 
as in figure|2] As an illustration, we also show (crosses) the BCGs 
in the SDSS/2MASS (Jci-matched subsample whose line of best fit 
(see below) is closest to the mean of the 1000 subsamples. 

There is a clear trend in each panel of figure |7] indicating 
that clusters with large velocity dispersions, and presumably large 
masses, tend to have brighter BCGs. A Spearman rank correlation 
coefficient analysis indicates that all these correlations are signifi- 
cant at a level greater than 99.9%. 

To explore whether the luminosity-aci correlation evolves 
with redshift, we analysed the intermediate-redshift EDisCS/ABK 
clusters and the low-redshift SDSS ones separately. We define L^ 
as the rest-frame ii'-band luminosity of the BCGs corrected ei- 
ther assuming no luminosity evolution (N.E.) or passive luminos- 
ity evolution for Zform = 2 or 5. If we parameterize the trends 
shown in figure |7] as L^ = Ca'^i, simple least-squares fits to the 
EDisCS/ABK data give 



N.E. 


a 


= 0.35 ±0.14 


orni ^ ^ 


a 


= 0.36 ±0.14 


orin ^ 


a 


= 0.34 ±0.14 



and for the mean of the 1000 SDSS/2MASS ad-matched cluster 
subsamples we get a = 0.35 ± 0.10 in all cases. 

Note that the correlations found and their estimated slopes are 
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Figure 5. Histograms showing tlie cluster velocity dispersion a^i distribu- 
tion of the model clusters. The solid line shows the distribution of the model 
CTcl values at z = 0.04. The dotted and dashed lines show the a^i distribu- 
tions at 2 = 0.36 and z = 0.76 respectively. The a^i values were selected 
using the second method described in the text with a^i = 200 km s^^. 



robust against changes in the way we account for the passive lumi- 
nosity evolution of the BCGs. 

The slopes of the Lk-CcI relations are the same for the low- 
and high-z samples. This seems to indicate that in the the range of 
redshifts explored here the rate of the BCG buildup and that of the 
cluster are linked in a way that does not depend strongly on cosmic 
time. 

Another interesting result is that, although significant, the de- 
pendency of I/^ (presumably proportional to the BCG's stellar 
mass A/bcg) o" '^ci is relatively weak. If we assume that the clus- 
ter mass Aid is proportional to R20o<j'ii, where i?200 is the radius 
of a sphere with interior mean density 200 times the critical density 
of the universe at the cluster redsh ift. If we also assume that, at a 
given redshift, _R200 oc CTcI (Carlbe rg et alJI 19971 : [Finn et alj|2005h . 
then A/ci oc crj?i, and M^qq (x m012±o.03 jj^^^,^ ^^^ stellar mass 
of BCGs changes only by ~ 70% over a two-order-of-magnitude 
range in cluster mass. This is not surprising: it has long been known 
that BCG luminosities measured inside fixed metric apertures pro- 
vide reaso nable standard ca ndles almost independent of the cluster 
richness ( Postman & Lauer 199 ^. 

Very recendv. lBrough et aj] ( l2008h and lStott et all J2008l) have 
found similar results for X-ray selected cluster samples. 



3 COMPARISON WITH GALAXY FORMATION 
MODELS 

In this section, we compare our emp i rical r esults with theoretical 
predictions from lOe Lucia & Blaizod 120071). These authors use a 
combination of large N-body coUisionless cosmological simula- 
tions and semi-analytic techniques. The semi-analytic models are 
grafted onto the hier archical merging tree s extracted from the Mil- 
lennium Simulation JSpringel et al.ll2005h . In this study, we com- 
pare our results with two versions of the semi-analytic model which 
differ for the prescriptions adopted for the supemovae feedback. In 



the following we will refer as model 1 to a model which employ 
the same supemovae feedback model used in lOe Lucia et alj 120061) 
where the ejection rates are computed on the basis of energy conser- 
vation arguments. We will refer as model 2 to a model which adopts 
the same supemovae feedback model as in ICroton et al.l 1200q) 
where the ejection r a tes are proportional to the star formation rate. 
lOe Lucia & BlaizotI 12007h have shown that model 1 results in a 
slower evolution of the K-band magnitudes as a function of red- 
shift, albeit producing K-band magnitudes of local BCGs brighter 
than a factor of about 0.8 mag with respect to results from model 
2. 

In order to compare our results to model predictions, we ex- 
tracted 90 haloes within the simulation box, uniformly distributed 
in log(mass) between 5 x 10^^ Mq and 5 x 10^^ Mq. We extracted 
the haloes from the output of the simulation at jz = and followed 
their history back in time by tracking, at each time, the main pro- 
genitor branch defined as in De Lucia & Blaizot ( 2007). For each 
SN feedback model, we have then created tables with rest-frame 
and observed-frame magnitudes in 5 different bands (B, V, R, I, 
J, and K) for the brightest galaxy in each halo. For each halo we 
have also calculated the 3D and projected velocity dispersions us- 
ing all galaxies within 2 x R200 at three redshifts (z = 0.04, 0.36, 
and 0.76). In order to simulate the errors associated with the ob- 
served CTci values, for each halo and each redshift we extracted a 
random number of objects, between 20 and 70 from all galaxies 
within 2 x R200 and computed the ID velocity dispersion and their 
confidence limits using the same method that was used to calcu- 
late th e velocity dispersions for the EDisCS clusters dHallidav et al.l 
[2OO4I1J. This method tries to mimic the observational uncertainties 
in the model predictions. 

When selecting the model cluster sample to compare with 
observations, there are two obvious alternatives. First, one could 
follow the evolution of the same model clusters, i.e., select clus- 
ters at all redshifts which, at a given redshift (e.g., z — 0), have 
o"ci > cTiim- Second, one could follow similar cluster "samples", 
i.e., select clusters that, at each redshift, have ad > crum- Obvi- 
ously, the first selection tells us what happens to the same set of 
galaxies, while the second one is probably closer to observations 
since it follows comparable samples. We will therefore follow the 
second selection method with criim — 200 km s^^ (64-70 model 
clusters, depending on redshift) and criim — 400 km s^^ (28-45 
model clusters, also depending on redshift), which seem reasonable 
limits for comparison with our dataset (cf. figure [TJ. Note that the 
(Tiim for EDisCS is not well defined since velocity dispersion was 
not explicitly part of the EDisCS cluster selection criteria. How- 
ever, since the results do not depend on the precise value of this 
limit, the choice is not critical. 

Figure [5] shows the model ad distribution at each redshift bin 
obtained using this selection technique. Although the model ad 
values were not explicitly matched in each redshift bin, as was done 
for the observational data, their distributions are well matched at all 
redshifts considered. Moreover, they are also well matched to the 
(Tci values of the observational data (cf. figure[T) 

The first result of the comparison is that the colours of the 
model galaxies are similar to the observed ones for model 1, i.e.. 



^ In the smallest haloes there are cases where the number of objects is 
less than 20. In these cases we take all objects for the calculation of the 
velocity dispersion. Note, however, that this situation seldom arises for the 
final comparison sample, once a minimum velocity dispersion cut is applied 
(as we do in our analysis). 
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Figure 6. Plot of log(M/Mo) versus log{l + z) for the semi-analytic 
models. Mq is the mean stellar mass of the BCGs at 2 = 0.04, and M the 
mean stellar mass of the BCGs at different redshifts. The triangles represent 
the data from model 1 with ciim = 200kms~^. The doted line is a least- 
squai'es fit to these data points. See legend for the meaning of the rest of the 
symbols and lines. 



close to those of passively-evolving stellar populations formed 
early. This suggests that even if the buildup of the galaxies' stel- 
lar mass is gradual, this growth is dominated by the accretion of 
old stellar populations rather than by new star formation. Model 2 
predicts somewhat bluer colours. 

Both models 1 and 2 predict a growth in stellar mass that is sig- 
nificantly stronger than observed: the stellar mass of model BCGs 
grows by about a factor of ~ 3-4 between z ~ 1 and z ~ 0, while 
the observations are compatible with no growth at all and seem 
to rule out a factor of two growth since z ~ 1 (cf. section 12.2b . 
The stellar mass growth in the models does not depend strongly on 
which feedback method is used, although it is marginally stronger 
in model 2. The stellar mass growth is also quite insensitive to the 
value of (Tiini, which agrees with the observed independence of the 
stellar mass growth on CTcI- These results are illustrated in figure[6] 

Finally, we compare the observed dependency of the BCGs' 
corrected K-haaA luminosity (or stellar mass) with the cluster ve- 
locity dispersion (cf. section [2J4}. For model 1, regardless of the 
exact value of (Tum, the models predict, at all the redshifts consid- 
ered, Ui- oc <Th " or M^c.a oc M] 



. Model 2 predicts L%^ oc a^^^ 
or MgQQ oc M^i''. The data suggest a much weaker dependency 
of the stellar mass on the cluster velocity dispersion or mass. That 
is clearly shown in figure [T] the dashed line represents the aver- 
age slope from the model predictions, which is significantly steeper 
than the data. 

One major caveat in the comparison of our observations with 
the semi-analytical model results is that we have used magnitudes 
measured inside a fixed metric aperture while the models compute 
total magnitudes. Observationally it is very difficult (virtually im- 
poss ible) to measure accu rate total magnitudes for the BCGs (see, 
e.g-. lGonzalez et alj2005h . Indeed, it is not even clear whether one 
should incl ude the intracluster l ight (ICL) as part of the BCG's 
luminosity JGonzalez et al.ll2007n . A 37kpc diameter aperture en- 



compasses only 25%-50% of the total light contained in the BCG 
and ICL (Gonzalez, Zaritsky and Zabludoff, private communica- 
tion). Interestingly, using the /-band data published by these au- 
thors we find that although their 37kpc aperture magnitudes show 
a dependency with Od compatible with our measurements, their to- 
tal BCG-l-ICL luminosities show a significantly steeper dependence 
(^^BCG+iCL ^ ^0.6-) Although the slope is still considerably shal- 
lower than the model predictions, this could explain at least part of 
the disagreement between models and observations. Note, however, 
that the scatter is large and their sample relatively small. 

On the other hand, the models do not have spatial informa- 
tion regarding the distribution of the BCG light, so aperture mag- 
nitudes cannot be calculated. It is also not clear whether the lu- 
minosities of the model BCGs contain some of what, observation- 
ally, we could call intracluster light. The luminosity of the model 
BCG includes the stars formed in the central cluster halo and the 
stars formed in galaxies that merged with the central galaxy, but not 
stars stripped from other cluster galaxies due to tidal and "harass- 
ment" effects, since these are not modeled. Therefore, any discrep- 
ancy (and/or agreement) found in the behaviour of the observed and 
model BCGs should be taken with caution, since not all the physi- 
cal effects which potentially contribute to the growth of the BCGs 
are taken into account. 

A possible way of reconciling our observations and the model 
predictions could be that, if there is a large stellar mass growth in 
BCGs (e.g., a factor of 4 since z ~ 1), these stars could lie out- 
side of the 37kpc apertures we have used in our measurements. 
In that case, such a growth would really represent the growth of 
the intracluster component rather than the central BCG. How such 
growth would depend on the cluster properties is difficult to pre- 
dict, but it is clear that the intracluster component could play a very 
important role in t he formation and evolution of clusters (see, e.g., 
iLin & Mohnl2004l) . For rec ent theoretical studies of the formation 
of the intracluster light, see lMonaco et al.l ( 120060 and lConrov et al.l 
( l2007h . 



4 CONCLUSIONS 

In this paper we present new iC-band data for 21 brightest 
cluster galaxies (BCGs) from the ESO Distant Cluster Survey 
(EDisCS). Combinin g these data with the photom etry published 
by lAragon-Salamanca. Baugh & KauffmannI ( Il998h we have built 
a sample of 42 BCGs with good quality near-infrared photometry 
in the 0.2 < 2 < 1 redshift range. We have also put together a low- 
re dshift {{z) = 0.07) compari son sample from the BCG catalogue 
of Ivon der Linden et al.l (12007) for which optical and near-infrared 
photometry is available from the SDSS and 2MASS. By compar- 
ing the properties of the BCGs in the low- and high-z datasets and 
those of their parent clusters we have carried out a detailed study of 
the evolution of these galaxies. Since, as we demonstrate, the prop- 
erties of the BCGs depend upon the velocity dispersion CTci (and 
thus presumably the mass) of the cluster they inhabit, we carefully 
match the ad distributions of our samples before comparing them. 
In this work the luminosities of the BCGs are measured inside a 
fixed metric circular aperture with 37 kpc diameter. The main con- 
clusions of our study are: 

(i) In agreement with previous studies, we find that the A'-band 
Hubble diagram for BCGs exhibits very low scatter (~ 0.35 mag) 
over a redshift range of < z < 1. 

(ii) The colour evolution of BCGs as a function of redshift 
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Figure 7. Plots of the residuals from the lines of best fit in figure |2] vs. logiQ(o-(.i). The top panel shows the residuals after subtracting the no evolution 
prediction. The middle and bottom panels show the same but for the passive evolution models with Zform = 2 and Zform = 5 respectively. The EDisCS 
clusters have been plotted as plus signs, the points from the ABK sample for which we have velocity dispersions are plotted as circles and the data for a 
typical SDSS/2MASS crd-matched subsample have been plotted as crosses. The SDSS/2MASS subsample shown has a slope very close to the mean slope of 
all 1000 subsamples (see text for details). The dotted Unes represent the least-squares fits to the EDisCS/ABK data and the solid lines represent the average 
least-squares fits to the 1000 SDSS/2MASS subsamples. The dashed lines represent the average slope predicted by the semi-analytic models described in the 
text. 



shows good agreement with population synthesis models of stellar 
populations that formed at z > 2 and evolved passively thereafter. 

(iii) In contrast with some previous studies (e.g., Aragon- 
Salamanca et al. 1998) we find that the rest-frame A'-band lumi- 
nosity evolution of the BCGs is also well reproduced by passive 
evolution models. The main reason for this discrepancy is the use 
of a different, hopefully more accurate, low redshift comparison 
sample with well-matched a^i distribution. 

(iv) The agreement of the observations with passive luminos- 
ity evolution models means that we do not detect any significant 
change in the stellar mass of the BCG since z ~ 1. Formally, we 
rule out a growth in the stellar mass larger than a factor of 1.9 be- 
tween z — 1 and z — (3-cr upper limit). This result does not seem 
to depend on the velocity dispersion of the parent cluster. 

(v) We find that there is a correlation between CTci (the ID ve- 
locity dispersion of the clusters) and the ii'-band luminosity of the 
BCGs (after correcting for passive evolution). Clusters with large 
velocity dispersions, and therefore masses, tend to have brighter 
BCGs, i.e., BCGs with larger stellar masses. This dependency, al- 
though significant, is relatively weak: the stellar mass of the BCGs 
changes only by ~ 70% over a two-order-of-magnitude range in 



cluster mass. Furthermore, this dependency doesn't change signifi- 
cantly with redshift. 



We have compared our observational results with the hi- 
erarchical galaxy fo rmation and evolution model predictions of 
lOe Lucia & BlaizotI J2007l) . One major difficulty in this compari- 
son is that we have measured magnitudes inside a fixed metric aper- 
ture while the models compute total luminosities. We find that the 
models predict colours which are in reasonable agreement with the 
observations because the growth in stellar mass is dominated by 
the accretion of old stars. However, the stellar mass in the model 
BCGs grows by a factor of ~ 3-4 since z = l,a. growth rate which 
seems to be ruled out by the observations. The models predict a 
dependency between the BCG's stellar mass and the velocity dis- 
persion (mass) of the parent cluster in the same sense as the data, 
but the dependency is significantly stronger than observed. It there- 
fore seems as if building realistic BCGs is another challenge for 
this kind of models, indicating that the underlying physics is far 
from being completely understood. 
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